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Abstract. 

Surface-plasmon-polariton (SPP) wave propagation guided by a metal slab in 
a periodically nonhomogeneous sculptured nematic thin film (SNTF) was studied 
theoretically. The morphologically significant planes of the SNTF on both sides of 
the metal slab could either be aligned or twisted with respect to each other. The 
canonical boundary-value problem was formulated, solved for SPP-wave propagation, 
and examined to determine the effect of slab thickness on the multiplicity and the 
spatial profiles of SPP waves. Decrease in slab thickness was found to result in more 
intense coupling of two metal/SNTF interfaces. But when the metal slab becomes 
thicker, the coupling between the two interfaces reduces and SPP waves localize to 
one of the two interfaces. The greater the coupling between the two metal/SNTF 
interfaces, the smaller is the phase speed. 
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1. Introduction 

Quasi-particles formed collectively by plasmons and polaritons bound to the interface of 
a metal and a dielectric material are called surface plasmon-polaritons (SPPs) [1, 2, 3, 4]. 
Although SPPs are quantum objects, trains of SPPs traveling along the interface are 
classically viewed as SPP waves. Much recent literature is devoted to the exploitation 
of SPP waves for diverse technological purposes, including sensing, imaging, and 
communications [5, 6, 7, 8]. 

This technoscientific ferment is due to a resonance phenomenon that arises when 
the energy carried by photons in the dielectric material is transferred to free electrons 
in the metal at that interface, and vice versa. Different dielectric materials will become 
differently polarized on interrogation by an electromagnetic field, thereby enabling a 
widely used technique for sensing chemicals and biochemicals [5, 7]. Furthermore, SPP 
imaging systems are used for high-throughput analysis of biomolecular interactions — 
for proteomics, drug discovery, and pathway elucidation [9, 10]. SPP-based imaging 
technology has been successfully applied to the screening of bioaffinity interactions 
with DNA, carbohydrates, peptides, phage display libraries, and proteins [11]. SPP- 
based imaging techniques are also going to be useful for lithography [2, 12]. Finally, 
as SPP waves can be excited in the terahertz and optical regimes, they may be useful 
for high-speed information on computer chips [13]. Whereas conventional wires are 
very attenuative at frequencies beyond a few tens of GHz, ohmic losses are minimal for 
plasmonic transmission [8] which enables long-range communications [14]. 

Even a cursory perusal of literature will show that only one SPP wave at a 
specific frequency is excited along the interface, if the partnering dielectric material 
is homogeneous. But if this material were to be periodically nonhomogeneous in the 
direction normal to the interface, the number of SPP waves — of different phase speeds, 
attenuation rates, and field structures, but the same frequency — can be higher. The 
material can be isotropic [15] or anisotropic [16, 17], the latter obviously providing 
more options due to birefringence [18]. The possibility of exciting multiple SPP waves 
provides exciting prospects for enhancing the scope of the applications of SPP waves. 
For sensing applications, the use of more than one distinct SPP waves would increase 
confidence in a reported measurement; also, more than one analyte could be sensed 
at the same time, thereby increasing the capabilities of multi-analyte sensors. For 
imaging applications, the simultaneous creation of two images may become possible. For 
plasmonic communications, the availability of multiple channels would make information 
transmission more reliable as well as enhance capacity. 

These possibilities are promised by the use of sculptured thin films (STFs) 
as the partnering dielectric material. STFs are anisotropic and unidirectionally 
nonhomogeneous thin films whose permittivity profile has been nanoengineered during 
the physical vapor deposition process [19, 20]. The possibility of multiple SPP waves 
localized to the single interface of a metal and a periodically nonhomogeneous sculptured 
nematic thin film (SNTF) has been demonstrated both theoretically [17, 21, 22] and 
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experimentally [23]. An SNTF is a special type of an STF — with a permittivity 
profile which is sculptured, during physical vapor deposition, only in one plane called 
the morphologically significant plane [20]. Also, SPP wave propagation along the 
interface of a metal and a chiral STF has also been shown both theoretically [24] and 
experimentally [25] to admit more than one type of SPP waves. 

A further increase in the number of SPP waves would require the use of multiple 
parallel metal/dielectric interfaces, which is already established well with isotropic 
dielectric materials [26, 27, 28]. Pursuing this line of thinking, we decided to solve 
the canonical problem of wave propagation localized to a sufficiently thin metallic slab 
inserted in an SNTF. As our supposition was confirmed in a preliminary study [29], we 
conducted a full-scale investigation, the details and the results of which are reported 
here. In Sec. 2, we present the theoretical formulation of the canonical boundary-value 
problem, where a dispersion equation is obtained. Representative numerical results are 
discussed in Sec. 3, and concluding remarks are given in Sec. 4. 

An exp(— iu>t) time- dependence is implicit, with u denoting the angular frequency, 
t the time, and i = y/—l. The free-space wavenumber, the free-space wavelength, and 
the intrinsic impedance of free space are denoted by k = u^/e fi , A = 2n/k Q , and 
Vo = \/ Ho/ e Q) respectively, with //,, and e being the permeability and permittivity of 
free space. Vectors are in boldface, dyadics are underlined twice, column vectors are in 
boldface and enclosed within square brackets, and matrixes are underlined twice and 
square-bracketed. The asterisk denotes the complex conjugate, and the Cartesian unit 
vectors are identified as u x , u y , and u 2 . 

2. Canonical Boundary- Value Problem 

Suppose that the region L_ < z < L + is occupied by an isotropic and homogeneous 
metal with complex- valued relative permittivity scalar e me t- The thickness of the metal 
slab is denoted by L met = L + — 

The regions z ^ L± are occupied by the chosen SNTF with periodically 
nonhomogeneous permittivity dyadic [17, 21] 

gs^)= e ^)%{z)-e^ f (z)-^\z)-^\^), z^L ± , (1) 

where the locally orthorhombic symmetry is expressed through the diagonal dyadic 

£ 0) = e a (z) u z u z + e b (z) u x u x + e c (z) u y u y (2) 

=ref 

and the local tilt dyadic 

g y {z) = (u x u x + u z u z ) cos \x{z)] + (u z u x - u x u z ) sin [x{z)\ + u^u^ (3) 
expresses nematicity. Both the relative permittivity scalars £ a ,b,c{z) and the tilt angle 
x{ z ) are supposed to have been nano-engineered by a periodic variation of the direction 
of the vapor flux during fabrication by physical vapor deposition [20, 23]. This periodic 
variation is captured by the vapor incidence angle [23] 

z^L ± , (4) 



Xv[Z) 



Xv ± $v 



Sill 



ti(z - L±) 
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that varies sinusoidally with z. The third dyadic in Eq. (1) was chosen as 



§JCT^) = + VLyXiy) cos7 ± + {xiyVL x - VL x VLy) sin7 ± + u z u z , (5) 

so that plane formed by the unit vectors u 2 and u x cos ^ + u y sin ^ is the 
morphologically significant plane for z ^ L±. Thus, there is sufficient flexibility in 
the formulation with respect to the twist 7 + — 7~ of the two morphologically significant 
planes. Note that 7 + = 7~ in our preliminary report [29], but that restriction was 
removed for the full-scale investigations reported here. 

Without loss of generality, let us choose the direction of SPP-wave propagation in 
the xy plane to be parallel to the x axis. Accordingly, we set 



E(r) = e(z) exp (inx) 
H(r) = h(z) exp (ikx) 



(6) 



where k is a complex-valued scalar. 

The axial field components e z (z) and h z (z) can be expressed in terms of the column 
vector 



[ f (»] = [e x (z) e y (z) h x (z) h y (z)\ 



(7) 



H 



via 



h z (z) 


where one 4x4 matrix 



A ± (z 



\A±(z)].[f(z)}, 
\A™\z)] ■ [f(z)] 



z^L±, 
z e (L_,L+) 











0^0 





U)£Q E a {z)e h {z) 







e d {z) [e a (z) -e b {z)\ . 

+ K r~{ — r~\ sm \x( z )\ cos [x{z)\ 

e a {z) e b {z) 



involves the auxiliary quantity 

£d{z) = £ a (z) e h {z)/ {e a (z) cos 2 [x(z)] + e b (z) sin 2 [x(z)}} 



cos7 ± 


sin7 ± 





























— sin7 ± 











cos7 ± 



and the other 4x4 matrix 



[A met (z)] 




























(8) 



(9) 



(10) 
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The column vector [f(z)\ satisfies the matrix differential equations 
d 



^ [f (*)] = i\PHz)] ■ [f(*)l , z^L ± 



and 



d_ 

dz 



[f(z)]=i[P met (z)]-[f(z)] , ze(L_,L_ 



(12) 
(13) 



where the 4x4 matrixes 

[E+(z)] = [AH*)] + 




e [e c (z) - e d {z)\ cos7 ± sin7 ± 
e [e c (z) sin 2 7 ± + e d (z) cos 2 7 ± ] 



to 






-£„ |£ c (z) COS 2 7 ± 



„ [e c (z) cos 2 7 ± + e rf (z) sin 2 7 ± ] 
e [e c (z) - e d (z)] cos7 ± sin7 ± 



fl 

-//„ o 







;i4) 



and 



\E met (z)] = [A met (z)] + 












/'•(, 








-A*o 




























(15) 



Equation (13) can be solved straightforwardly to yield 
[f (L + )] = exp {z \P met ] (L+ - L_)} • [f(L_] . 



(16) 



Equation (12) requires numerical solution by the piecewise uniform approximation 



technique [22, 30]. The ultimate aim is to determine the matrixes Q^ 1 
the relations 



[f (L± ± 2ft)] 



that appear in 
(17) 



to characterize the optical response of one period of the SNTF on either side of the 
metal slab. Basically, this technique consists of subdividing each period of the SNTF 
into a cascade of electrically thin sublayers parallel to the plane z — 0, and assuming 
the dielectric properties to be spatially uniform in each sublayer. A sufficiently large 
number N + 1 points z^ = L± ± 2ft (n/N), n G [0, N], are defined on each side of the 
metal slab and the matrixes 

2ft 



exp < ±i 



Z n-1 + Z n ^ 



N 



n e [1,N] 



18) 



are calculated for a specific value of k; then 

gi = ogi • [mi J •••••[£]• tEf] • (is) 

A sublayer thickness 2Q/N = 2 nm was adequate for the results reported in Sec. 3. 

By virtue of the Floquet-Lyapunov theorem [31], we can define the matrixes [Q 
such that 



[Q±] = exp { ±i2tt[Q ± ] 



(20) 
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Both [Q^] and [Q^\ share the same eigenvectors, and their eigenvalues are also related 
as follows. Let [t*] , n G [1,4], be the eigenvector corresponding to the nth eigenvalue 



cr± of 



then, the corresponding eigenvalue of [Q~"~] is given by 
Mat 



Q- 



2Q 



(21) 



The electromagnetic fields of the SPP wave must diminish in magnitude as z — > 
±oo. Therefore, in the half-space z > L + , we first label the eigenvalues of [Q ] such 
that Im[a^ 2 ] > and then set 

~ Af 



lHL + )]= [t 



(1) [ t +](2) 



At 



(22) 



where Af 2 are unknown scalars; the other two eigenvalues of [Q } describe fields that 
amplify as z — > +oo and cannot therefore contribute to the SPP wave. A similar 
argument for the half-space z < L_ requires us to ensure that Im [aj" 2 ] < and then to 
set 



lHL-)]= [t 



[t 



-1(2) 



A: 



(23) 



where A± 2 are unknown scalars. 

Combining Eqs. (16), (22), and (23), we obtain a matrix equation which may be 
rearranged as 



[M(«)l 



'Af 




" " 


At 







A~, 







. A 2 . 








(24) 



For a nontrivial solution, the 4x4 matrix [M(W)] must be singular, so that 

det [M(«)] = (25) 
is the dispersion equation for SPP-wave propagation. 



3. Numerical Results and Discussion 

A Mathematica™ program was written and implemented to solve (25) to obtain k 
for specific values of 7 + and 7". The dispersion equation (25) was solved using the 
Newton-Raphson technique [33] for three different values of the twist between the 
morphologically significant planes on either side of the metal slab; we chose 

(i) 7 - = 7 +, 

(ii) 7 - = 7 + + 90°, and 

(iii) 7- = 7 + + 45°, 
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while 7 + was kept as a variable. For each choice, the boundaries of the metal slab 
were taken to be at L± = ±7.5, ±12.5, ±25, or ±45 nm. These selections adequately 
represent the results of our investigation. 

The free-space wavelength was fixed at A = 633 nm for all calculations. The metal 
was taken to be bulk aluminum [e me t = —56 ± 21i). The skin depth of aluminum at 
the chosen wavelength is A met = {im [/c ov /e met ] } = 13.24 nm, a quantity of interest 
in relation to the thickness of the metal slab. The SNTF was chosen to be made of 
titanium oxide [29, 32], with 

e a (z) = [1.0443 ± 2.7394u(z) - 1.3697w 2 (^)] 2 ' 

e b (z) = [1.6765 + 1.5649v(z) - 0.7825w 2 (z)l 2 

} , (26) 

e c (z) = [1.3586 + 2.1109v(z) - 1.0554i; 2 (z)] 2 

X(z) = tan" 1 [2.8818 tan x^)] 

where v(z) = 2x v {z)/n. The angles Xv and 5 V were taken to be 45° and 30°, respectively, 
for all results presented here. 



2.05 
2.00 



' ■ ■ 1 


:■■ ■ 


■ ■ 


■ ■ 




































































-«-«• 


-■ 























y + (deg) 



(a) 



■ ■ 1 






-m-m- 


i ■ ■ i 


I ■ ■ 


• • 




♦ * 


< • • 


• •- 
















lA* 1 

■-▼ T 
























■ 


1 


--- 

















I** 


■ ■ i 






























— 














-i 














tA-A j 

< 












15 30 45 60 75 90 

y + (deg) 


■ ■ 


■ 


m-m 


mm- 


i ■ ■ 


i ■ ■ i 
















Mi 






















.ill 


k-A A 


-■ 






























' ▼ ▼ ' 



0.025 
0.020^P 
0.015 £ 



30 45 60 75 90 15 30 45 60 75 

Y + (deg) (M y + (deg) 



1.95 
1.90 





■ ■ 


■ ■ i 


i ■ ■ i 


1 ■ ■ 1 


















































• • 1 








A A j 














































i A Aj 


,A A< 


rill 





























0.0075'-=' 

E 



15 30 45 60 75 90 15 30 45 60 75 90 

y + (deg) ( C ) y + (deg) 



0.10 


2.50 


09 
0.020 


2.40 


0.015 ^ 

E 


2.10 

""6 

^ 2.05 
01 


0.010 


CH 2.00 

1.95 


0.005 


1.90 


0.000 


1.85 



1 


1 ■ ■ 




■-■ 


■ ■ I 


■ ■ i 




















































4 









15 30 45 60 75 90 

T + (deg) 



(d) 



15 30 45 60 75 90 

y + (deg) 



Figure 1. Variation of real and imaginary parts of n/k with 7 + , when 7 = 7+. (a) 
L ± = ±7.5 nm, (b) L± = ±12.5 nm, (c) L± = ±25 nm, and (d) L± = ±45 nm. 
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Figure 2. Variation of the Cartesian components of the time-averaged Poynting vector 
P(x,z) along the z axis when x = 0, L± = ±7.5 nm, and 7~ = 7 + . (a-c) 7 + = 0, 
and (d-f) 7+ = 25°. (a) n/k = 2.6387 + ^0.1839, (b) n/k = 2.0964 + ^0.009997, (c) 
n/k = 1.9048 + i0.02696, (d) n/k = 2.6399 + i0.1848, (e) n/k = 2.09285 + ^0.00988, 
and (f) n/k Q = 1.9103 + iO. 02405. The x-, y- and z-directed components of P(x, z) arc 
represented by solid red, dashed blue, and chain-dashed black lines, respectively. 



3.1. 7- = 7+ 

Let us begin with the case when the morphologically significant planes of the SNTF on 
both sides of the metal slab are aligned with each other and make an angle 7" = 7" 1 " 
with respect to the direction of SPP-wave propagation in the xy plane. This case had 
been briefly studied in a predecessor communication [29] , but is now discussed in detail. 

The real and imaginary parts of n/k which satisfies Eq. (25) are presented in Fig. 1 
as functions of 7 + G [0°,90°]; by virtue of symmetry, the solutions for 180° + 7 + and 
360° — 7 + are the same as for 7 + . For the thinnest metal slab (L± = ±7.5 nm), the 
solutions are organized in five branches which span the entire range of 7 + . As the metal 
slab thickens to 25 nm (L± = ±12.5 nm), the number of branches does not change, 
but only four of those branches span the entire range of 7" 1 " and one is confined to 
7 + G [0°,49°]. Further thickening of the metal slab to 50 nm (L± = ±25 nm) leads to 
five values of k satisfying the dispersion equation only in the range 7 + G [0°,35°], four 
in the range 7 + G (35°, 37°], and three in the range 7 + G (37°, 90°]. Finally, for a 90-nm 
thick metal slab (L± = ±45 nm) only three solutions exist for 7 G [0°, 36°] and two for 
7 G (36°, 90°], these solutions being the same as for a single metal/SNTF interface [22]. 
We conclude that, as the thickness of the metal slab is increased, the coupling between 




Figure 3. Same as Fig. 2 except for L± = ±45 nm. (a) K,/k = 2.4549 + i0.04173, (b) 
«/fc = 2.08034 + ^0.003574, (c) n/k = 1.8683 + ^0.00734, (d) n/k = 2.4557+i0.04181, 
(e) n/k = 2.0773 + ^0.00363, and (f) n/k = 1.8830 + ^0.00413. 



two metal/SNTF interfaces z = L_ and z = L + decreases; ultimately, the two interfaces 
decouple from each other when the thickness L met significantly exceeds twice the skin 
depth A met in the metal. 

The solutions in Fig. 1 can be categorized into three sets. The first set comprises 
those solutions for which Re[«/fc ] lies between 2.3 and 2.7. This set has two branches 
when the metal slab is 15-nm thick, both branches spanning the entire range of 7 + . As 
the thickness L met increases, the two branches come closer to each other and eventually 
merge completely. The second set comprises solutions for which Re lies between 
2.05 and 2.1. It also has two branches when L met = 15 nm, both branches coalescing 
into one branch as the thickness increases. Complete merger of the two branches of 
the first set occurs at a value of L met smaller than that for the two branches of the 
first set. Regardless of the value of L met , solutions in the second set can be found 
over the entire range of 7 + . The third set consists of solutions lying on just one branch 
(1.85 < Re[/c/A; ] < 1.95), but the maximum value of 7 + on this branch decreases rapidly 
from 90° as L met increases from 15 nm. 

The foregoing categorization is also meaningful as the spatial profiles of the fields 
are very similar for both solutions (for a specific value of 7 + ) in the first two sets. The 
same remark can be made for the spatial profiles of the time-averaged Poynting vector 

P{x,z) = ^Re[E(x,z) x H*{x,z)} . (27) 
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Representative plots of P(0, z) against z are presented in Figs. 2 and 3 for combinations 
of (i) three values of k, one from each set, and (ii) two values of 7 + . Whereas 7 + = 
was chosen because it represents the propagation of SPP waves in the morphologically 
significant plane of the SNTF on either side of the metal slab, 7 + = 25° was chosen for 
the direction of SPP-wave propagation oblique to that plane. 

Figure 2 shows the spatial profiles of P(0, z) for L± = ±7.5 nm. The chosen values 
of n/k Q are: (a) 2.6387 + z0.1839, (b) 2.0964 + z0.009997, and (c) 1.9048 + z0.02696 for 
7+ = 0; and (d) 2.6399 + z'0.1848, (e) 2.09285 + i0.00988, and (f) 1.9103 + i0.02405 
for 7 + = 25°. In order to determine these profiles, first Af was set equal to unity 
and then the remaining coefficients were found using Eq. (24) for all n/k — except for 
n/k = 2.0964 + iO. 009997, for which A% = 1 was set and the remaining coefficients 
were found using Eq. (24). The spatial profiles shown allow us to conclude that 

(i) the SPP waves are bound strongly to both metal/SNTF interfaces, and 

(ii) the power density mostly resides in the SNTF. 

Representative calculated values of the penetration depth A^ of the SPP wave 
into the metal, defined as the distance along the z axis (in the metal slab) at which 
the amplitude of the electric or magnetic field decays to e -1 of its value at the nearest 
metal/SNTF interface z = L±, are tabulated in Table 1 for L± = ±7.5 nm. The 
penetration depths for all SPP waves are ~ 12.7 nm, which confirms strong coupling of 
the two metal/SNTF interfaces. This is not surprising since A~ = A+ is very close to 
A me j . 

Table 1. Penetration depths A+ = A~ for L± = ±7.5 nm and -f~ = 7 + . The 
solutions are numbered in descending values of Rc[K/fc ]. 



A+ = A; (nm) 


7 + = 7 


Solution 


1st 


2nd 


3rd 


4th 


5th 


0° 


12.5458 


12.6558 


12.7780 


12.7882 


12.8564 


15° 


12.5456 


12.6557 


12.7786 


12.7886 


12.8555 


30° 


12.5450 


12.6553 


12.7801 


12.7896 


12.8532 


45° 


12.5443 


12.6548 


12.7821 


12.7907 


12.8500 


60° 


12.5436 


12.6543 


12.7839 


12.7916 


12.8470 


75° 


12.5431 


12.6539 


12.7852 


12.7921 


12.8454 


90° 


12.5430 


12.6538 


12.7857 


12.7922 


12.8451 



The spatial profiles of P(0, z) for L± = ±45 nm are presented in Fig. 3 for six 
selected values of n/k identified in the figure caption. These profiles were obtained after 
A\ was set equal to unity, except that A\ = 1 was used for n/k = 2.08034 + i0. 003574. 
The spatial profiles of P(0, z) for L± = ±12.5 nm and ±25 nm have not been shown 
here because the spatial profiles for these cases are similar to those shown in Figs. 2 and 
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3. Representative values of the penetration depths A^ for L± = ±45 nm are given in 
Table 2. The penetration depths for this case are of the same order as for L± = ±7.5 nm. 
As the thickness of the metal slab is 90 nm, it can be safely concluded that the SPP 
waves propagating on the two metal/SNTF interfaces are not coupled to each other but 
are propagating independently. 

Table 2. Penetration depths A+ = A~ for L± — ±45 nm and 7~ = j + . The solutions 
are numbered in descending values of ~Re[K/k ]. 



A+ = A~ (nm) 


7 + = 7 


Solution 


1st 


2nd 


3rd 


0° 


12.6206 


12.7843 


12.8691 


15° 


12.6205 


12.7848 


12.8669 


30° 


12.6201 


12.7860 


12.8608 


45° 


12.6196 


12.7877 




60° 


12.6192 


12.7891 




75° 


12.6188 


12.7901 




90° 


12.6187 


12.7904 
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Figure 5. Variation of the Cartesian components of the time-averaged Poynting 
vector P(x,z) along the z axis when x = 0, L± = ±7.5 nm, and 7~ = 7 + + 90°. 
(a-c) 7 + = 0, and (d-f) 7 + = 25°. The following values of n were chosen for rough 
correspondence with those in Fig. 2: (a) n/k Q = 2.3753 + iO. 005699, (b) n/k = 
2.09013 + i0.009135, (c) n/k = 1.9133 + ^0.004397, (d) n/k = 2.3756 + ^0.005713, (e) 
«/fc = 2.08797 + iO.009324, and (f) n/k = 1.9165 + i0.01117. 



3.2. 7" = 7+ + 90° 

The real and imaginary parts of n/k Q that satisfies the dispersion equation (25) for 
7^ = 7 + + 90° are given in Fig. 4 as functions of 7 + G [0°, 90°]. Due to the symmetry 
of the problem, the solutions for 90° ± 7+, 180° ± 7+, 270° ± 7+ and 360° - 7+ are 
the same as for 7 + . We found five solutions which span the whole range of 7" 1 " for 
L± = ±7.5 and ± 12.5 nm. However, when L± = ±25 and ±45 nm, only four solutions 
exist for 7 £ (37°, 53°), but five for other values of 7+ £ [0°, 37°] U [53°, 90°]. 

For the thickest slab considered, we see from comparing the results presented in 
Fig. 4(d) with those in Fig. 1(d) — which are same as that for a single metal/SNTF 
interface — that the two metal/SNTF interfaces are actually uncoupled from each other. 
The solutions in Fig. 4(d) represent SPP-wave propagation guided by a metal/SNTF 
interface with the direction of propagation in the xy plane either making an angle 7 + 
or 7~with the morphologically significant plane of the SNTF. 

As in Sec. 3.1, the solutions of the dispersion equation can be categorized into three 
sets with the same criterions as given in Sec. 3.1. However, in this case, the two branches 
in either of the first two sets do not merge into one branch as L met is increased. Instead, 




Figure 6. Variation of the Cartesian components of the time-averaged Poynting 
vector P(x, z) along the z axis when x = 0, L± = ±45 nm, and 7~ = 7 + + 90°. 
The following values of k and 7 + were chosen to highlight the uncoupling of the two 
metal/SNTF interfaces, when the metal slab is sufficiently thick, (a-e) 7 + = 25° and 
(f) 7+ = 65°. (a) n/k = 2.4558 + £0.04214, (b) n/k = 2.07725 + £0.003595, (c) n/k„ = 
1.8830 + £0.004085, (d) K /k = 2.45833 + £0.0424391, (e) n/k = 2.06775 + £0.00381762, 
and (f) n/k = 1.8830 + £0.004085. 



the two branches remain distinct, each holding for SPP-wave propagation guided by one 
of the two metal/SNTF interfaces uncoupled from the other metal/SNTF interface. The 
single branch in the third set actually vanishes for mid-range values of 7 + , with the two 
parts of that branch signifying independent propagation guided by the two metal/SNTF 
interfaces. 

Examination of the field profiles confirms the conclusions made in Sec. 3.1 regarding 
the effect of the thickness of the metal slab. The asymmetry in the alignment of the 
morphologically significant planes for z < L_ and z > L + is reflected in the spatial 
profiles of the time- averaged Poynting vector presented in Figs. 5 and 6 for L met = 15 
and 90 nm, respectively. Fig. 5 also suggests that the coupling of the SPP-wave is 
directly proportional to the real part of wavenumber. So slower SPP waves are more 
strongly coupled to both interfaces. This observation is also confirmed by the spatial 
profiles of SPP waves when L met = 25 and 50 nm (not shown). Fig. 6 shows the 
uncoupling of the two metal/SNTF interfaces when the metal slab is sufficiently thick. 
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Figure 7. Same as Fig. 1 except 7 = 7+ + 45°. 



3.3. 7 - = 7 + + 45° 

The solutions of the dispersion equation (25) for 7" = 7" 1 " + 45° are given in Fig. 7 for 
7 G [0°, 180°]. By virtue of symmetry, the solutions for 7 + + 180° are the same as for 
7 + . We found that five solutions exist for the entire range of 7 + for L± = ±7.5 nm. 
When L± = ±12.5 nm, five solutions exist for 7 G [0°,37°] U [98°, 180°] but only four 
for 7 G (37°, 98°). Further thickening of the metal to 50 nm (L± = ±25 nm) yields five 
solutions for 7 G [0°, 36°] U [99°, 180°] but four for 7 G (36°, 99°). Decoupling of the two 
metal/SNTF interfaces becomes very pronounced for L± = ±45 nm, when all solutions 
found are the same for either (i) a metal/SNTF interface for which the direction of 
propagation in the xy plane makes an angel 7 + with the morphologically significant 
plane, or (ii) a metal/SNTF interface for which the direction of propagation in the xy 
plane is inclined at 7" 1 " + 45° to the morphologically significant plane. 

Similar to Sees. 3.1 and 3.2, the solutions can be grouped into three sets, with the 
same criterions given in Sec. 3.1. Representative field profiles for 15-nm- and 90-nm- 
thick metal slabs are given in Figs. 8 and 9, respectively. We selected one value of k 
from each set at 7 + = 25° and 150°. Whereas 7 + = 25° corresponds to the propagation 
in the xy plane at an angle of 25° with respect to the morphologically significant plane 
in the region z > L + and at 70° in the region z < L_, the analogous angles are 30° 
and 15°, respectively, when 7 + = 150°. Fig. 8 shows that SPP waves are guided by 
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Figure 8. Variation of the Cartesian components of P(x, z) along the z axis 
when x = 0, L± = ±7.5 nm, and 7^ = 7+ + 45°. (a-c) 7 + = 25, and (d- 
f) 7+ = 150°. (a) n/k = 2.6423 + i0.1865, (b) n/k = 2.08764 + i0.009246, (c) 
re/fc = 1.9159 + ^0.009486, (d) n/k = 2.6398 + ^0.1847, (e) n/k = 2.09389 + i0.01000, 
and (f) K/k = 1.9108 + ^0.02442. 



both interfaces of thin metal slab, but the power density profile is asymmetric due to 
the fact that morphologically significant planes are not parallel to each other in the two 
regions occupied by the SNTF. Figure 9 shows that when L met = 90 nm, any SPP wave 
propagates predominantly guided by one of the two metal/SNTF interfaces. We deduce 
from these two figures that the conclusions drawn in Sees. 3.1 and 3.2 hold true for this 
case as well, and therefore are general enough. 

4. Concluding Remarks 

We formulated and solved a canonical boundary-value problem to examine the 
characteristics of surface-plasmon-polariton waves guided by a thin metal slab inserted 
in a periodically nonhomogeneous sculptured nematic thin film. The morphologically 
significant planes on the two sides of the metal slab could be either parallel to or twisted 
with respect to each other. When the metal slab is very thin, its two interfaces with the 
SNTF couple to each other, thereby generating more modes of SPP-wave propagation. 
As the metal slab thickness, the two interfaces decouple. Finally, the two interfaces of 
a sufficiently thick metal slab independently guide SPP waves. 
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Figure 9. Same as Fig. 8 except that L± = ±45 nm. (a) n/k Q = 2.4586 + ^0.04246, (b) 
n/k = 2.07725+^0.003595, (c) n/k = 1.8891+^0.002576, (d) n/k = 2.4559+^0.04227, 
(c) n/k = 2.07916 + ^0.003560, and (f) n/k = 1.8737 + i0.006143. 
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